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■^^ ■ Quarks in the adjoint representation have been a subject of study 

iq \ for both conceptual and practical purposes. Conceptually, their dif- 

ferences when it comes to confining and chiral symmetry properties 
has long been suspected to hold important information on the rela- 
tion of these two distinguished properties of QCD-like gauge theories. 
Practically, they have been studied as both a possibility to access fi- 
nite density quark systems as well as candidate theories for technicolor 
in beyond-the-standard-model settings. The most elementary object 
c^ ■ describing such particles is their propagator, though it being gauge- 

dependent. Its properties in the minimal Landau gauge are investi- 
gated here both in the quenched and unquenched case for a range of 
^. | lattice parameters using the Wilson formulation for the gauge group 

(■*") i SU(2). It is found that the propagator shows pronounced differences 

OQ | to the case of fundamental quarks, especially towards the chiral limit. 

■^t" ! 1 Introduction 

O 

QCD belongs to a wider class of theories, in which gauge fields of a given 
Lie algebra interact with quark fields in some representation of the corre- 
sponding Lie group. Out of these possibilities, it is quarks in the adjoint 
representation which are quite interesting for several reasons. 
C^ | In QCD it is found that in the quenched case the chiral and decon- 



finement transition coincides [l|. Even in the unquenched case, where the 
transition becomes a cross-over, both transitions remain at the very least 
close by [2H4J]. This has been found for all gauge groups investigated so 
far [5rl8(]. This fact has been intriguing since a long time, and there were 
a multitude of investigations to understand whether there is a fundamental 
relation between both effects. 

However, for adjoint quarks, i. e. quarks in the same representation as the 
gauge group, the situation is drastically different. There, both temperatures 
differ already in the quenched case by a gauge-group-dependent factor of four 



or more [9]. This remains true in the unquenched case 10J, [ll]. Again, this 
is true for all gauge groups investigated so far. The question is thus whether 
the relation of confinement and chiral symmetry breaking of the fundamental 
case is just purely coincidental, and not true for general representations, or 
whether adjoint quarks just work differently than fundamental quarks. That 
is particularly interesting as adjoint quarks do not have a sign problem at 
finite baryon density, and thus are an interesting proxy for QCD [12J, [l3| ■ 
However, this makes only sense, if the mechanisms involved are not (too) 
different. 

Once the adjoint quarks are dynamical, the theory shows again a very 
different behavior depending on the number of quark flavors, at least in the 
case of SU(iV) gauge groups. For a single Weyl fermion, the theory becomes 
supersymmetric. On the other hand, already at two flavors, the theory 



becomes likely quasi-conformal, i. e. walking, in the chiral limi t |14l. Il5l| . It 
is therefore highly interesting for technicolor phenomenology |l4-ll7l] , and 
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therefore has been investigated repeatedly 

To understand these differences, a logical starting point is the most ba- 
sic object describing these quarks: The quark propagator. For QCD, this 
propagator has been calculated many times with different methods, see e. 
g. [32 42]. It shows clear signs of dynamical chiral symmetry breaking, and 
there are hints that it also shows that the quark is an unphysical particle 
due to the absence of a Kallen-Lehmann representation, i. e. by positivity 



violation 431. However, so far there are very few investigations for adjoint 
quarks [la . 



Mm 

Here, this situation will be improved on, using lattice simulations, for 
the case of the SU(2) gauge group. Since the quark propagator is gauge- 
dependent, this has to be done in a fixed gauge, which will be the minimal 
Landau gauge [47] . After a brief discussion of the simulation details in sec- 
tion [21 first the quenched and then the unquenched case will be investigated 
in sections [3] and HI respectively. The properties of the quark propagator 
will be detailed for a range of mass values and lattice parameters, in both 
momentum and position space. This results will be summarized together in 
the closing discussion in section [5j These investigations continue and com- 
plement the analysis of the gluonic correlations functions of the unquenched 
theory in [3(| 

2 Lattice setup 

2.1 Configurations 

In the following the adjoint quark propagator will be investigated in the 
quenched and unquenched case. Since the unquenched configurations were 



obtained from the authors of [la . Il9l . l2ll . |48] , the quenched simulations were 



chosen to match the unquenched ones, as has already been done in [3 



Table 1: The quenched configurations employed. Nt is the temporal and 
N s the spatial extent of the lattice, and the total volume in lattice units is 
N t x Ng. The lattice spacing has been determined using the data of [491 ]. 
setting the string tension to (440 MeV) 2 . In all cases the quenched quark 
mass in physical units were m =0.01, 0.1, 0.5, 1, 2, and 10 GeV. 



p 


a" 1 [GeV] 


Lattice sizes 


Configurations 


2.221 


0.985 


8 4 , 16 4 , 24 4 


1035, 667, 293 


2.328 


1.31 


8 4 , 16 4 , 24 4 


828, 850, 225 


2.457 


1.97 


8 4 , 16 4 , 24 4 


1035, 850, 199 


2.656 


3.94 


16 4 , 24 4 


750, 138 



Since the unquenched ones are not made by us, we refer for the details of 
their creation to the original literature [48J], and will detail here only their 
matching quenched case, as well as the additional calculations of the quark 
propagator. 

Thus, in the quenched case the action was chosen to be the standard 
Wilson action, and configurations were created using the methods described 
in [50J, i. e. with a combination of heat-bath and overrelaxation sweeps. The 
list of the lattice and configuration creation parameters for the quenched 
configurations are given in table [TJ For the unquenched case, the list of 
configurations can be found in [30J. Here, all but the configurations with 
a time extent of Nt = 64 have been included, the latter being beyond the 
computational resources available to this project. For the quenched results, 
the usual QCD scale setting will be used 49 J . 50J|. However, for the dynamical 
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30j, |48|] . Therefore results will 



case, this is a non-trivial problem 18|, Il9l . 
be presented both in lattice units and in physical units chosen such that the 
lightest gau! ge-invariant ++ state, which appears to be the glueball [191 ] . 
will be given the mass of the, very likely, observed 5l|, ]52| Higgs boson, i. e. 
125 GeV, following the arg umentation on the association of the Higgs mass 
with physical states of 53j, |54]. This will be detailed more below in section 

m 

Since the quark propagator is gauge-dependent, it is necessary to fix a 
gauge. To be compatible with previous investigations 30j, |44|, and because 
of the general advantages 47], Landau gauge will be chosen. To cope with 
Gribov copies, this is further specified to the minimal Landau gauge [47J. 
However, at least for scalar particles no significant influence on the propaga- 
tor due to treatment of Gribov copies has been observed so far [55J], and thus 
this choice is expected to be of little relevance here. The actual gauge-fixing 
is then performed using the method described in [50]. 



2.2 The adjoint quark propagator and lattice corrections 

2.2.1 Continuum 

Due to the unbroken color symmetry the quark propagator in Landau gauge 
in the continuum in four Euclidean dimensions is parametrized by two scalar 
function A and B as 

cab ( s _ rab A (p 2 h^ + B(p 2 ) _ b 2 , 7/JV + M{p 2 ) 

{P) A(p 2 ) 2 p 2 + B(p 2 ) 2 {P ' p 2 + M{p) 2 {) 



Z(P 2 ) 
M(p 2 ) 



1 



A{ P y 

B(p 2 ) 
A(p 2 ) 



with the Euclidean Dirac matrices 7«. The wave- function renormalization 
Z and the mass function M represent an alternative parametrization of 
the dressing functions A and B. Especially, if a solution to the equation 
M(—m 2 ) = —m 2 exists the mass function determines the pole mass. How- 
ever, even if a solution to the equation exists, it cannot be directly accessed 
with just the knowledge of the domain p 2 > 0, except if it would hold that 
M(p 2 ) = M(—p 2 ), what is generically not the case. 

Instead, to identify a possible mass, the Schwinger functions 43u47|. 



1 l'°° Z(r> 2 ) 

A " = kJ dpC ° s{tp) p 2 + M(p 2 ) 2 (2) 

1 f°° , , ,Z(p 2 )M(p 2 ) 
^ = -jf dpco^ /^i , (3) 

where v refers to the vector part and s to the scalar part, can be used. 
Just as ordinary space correlators [56(, these functions will be exponentially 
decaying at large times, with the decay constant being the pole mass, if such 



a pole mass exists 47|. This will be investigated in sections 13.21 and 14.31 for 
the quenched and unquenched case, respectively. 

However, before this can be done, the propagator has to be renormal- 
ized, as it is logarithmically divergent in four dimensions. To renormalize 
the quark propagator requires in general two conditions, one for each ten- 
sor component, i. e. A and B. These will be chosen to be Ar([i)=1 and 
Br(h) = m, where m is the tree-level mass. In the continuum, this can 
be achieved by two multiplicative renormalization factor Za and Zb, which 
obey in Landau gauge in this scheme Za = Zb [33j], such that the M func- 
tion is renormalization-group-invariant. On the lattice, the situation will be 
different due to lattice artifacts. 



2.2.2 Lattice implementation 



Following [48] for the lattice Dirac operator the Wilson operator with anti- 
periodic boundary conditions in time, but periodic boundary conditions in 
space, will be used, 

W(x, y) = - l - J](l - ^)U^6 X+Iiy + (m + 4)5 xy . (4) 

U? are the links in the adjoint representation, which are obtained from the 
links U^ in the fundamental representation by 



^ = ^^V[/J, (5) 



with the Pauli matrices a a . The corresponding lattice quark propagator in 
momentum space is then given by 



^fn^J-ST fl —„( 2 -* P A\ ^-U\.T „, oSt, (I5E1 



D L = 5 ab 



Bl(p 2 ) [m + £„(l " cos l^f - A L (p*)iZ,J^m 



N u 



B L (P 2 ) 2 lm+ £ u (l " cos ( 2 -^)Y + A L (P 2 ) 2 (E u sin ^''" ^ ^ ' 



where P^ are the lattice momenta, which take the half-integer values 1/2, ..., 
(N + l)/2 in time direction. This is the only direction which will be con- 
sidered in the following, and therefore always all momenta will be tacitly 
taken to be in the time direction. The reasons for this is to make full use 
of the extended time direction of the unquenched configurations. In the 
quenched case, all lattices are symmetric, but the anti-periodic direction 
will still make a difference in a finite volume. To obtain Dl requires to 
invert and Fourier-transform the Wilson operator (jj]). This has been ac- 
complished in a standard manner using a bi-conjugate gradient inversion on 
a plane-wave source for each momentum point, and subsequent projection 
onto the conjugate plane-wave. This code has been thoroughly checked for 
the free case, and for the case of pure global color symmetry, before going 
to the interacting case. 

The ultimate aim is to determine an approximation to the continuum 
functions A and B. Thus, the following functions will be determined, where 
the momentum is already chosen to be in the time direction, 

tr7 £>L 
pA L -- 



B L 



(tr 70J D L ) 2 + {txD L Y 
trD L 



where p is the relevant physical momentum component, 2sm(7rP/Nt). Of 
course, in the continuum limit this reduces to a trivial identity, Al = A and 



Bl = B. On a finite lattice, both are affected by lattice artifacts. Already 



for the free case, the exact results are 56f] 



in. _ :„ f 27rP o \ 



Af cc = sin 



\ L J 



ofree , i /27tP 

B L 1CC = m + 1 - COS I — - — 

Thus, there is an additive, momentum-dependent lattice correction to the 
B function and a multiplicative one to the A function. The A function can 
be corrected by division by an appropriate factor. The B function can be 
corrected by either subtracting or normalizing the Bl function. The latter 
option will be implemented here, as the B function can then immediately 
be interpreted as the deviation from the free case. The corrected dressing 
functions are therefore 



.4 



A, 



/(free 



nfree ' *> ' 

respectively. 

The calculation of the Schwinger functions on the lattice afterwards is a 
straightforward discretization of the inverse Fourier transformations (J2H31), 
taking into account the whole momentum range [471 ]. 

Statistical errors are calculated using bootstrap with 1000 resamplings 
and permitting asymmetric errors [50j. For derived quantities, including the 
Schwinger functions, errors are propagated. 

3 Quenched 

For the quenched results, to comply with the QCD orientation of this inves- 
tigation, the scale is set by the string tension. It is given the value of (440 



MeV) , using the results of [49(] as input. 



3.1 Results in momentum space 

To appreciate the importance of the lattice corrections, the uncorrected 
functions Al and Bl are shown for several example lattice settings in figure 
[H as functions of the lattice momenta. Interestingly, the results do not 
deviate strongly from each other. In comparison, when including the lattice 
corrections, the situation changes, see figure [2j It is immediately visible 
that the independence of the mass on the regulator is not obtained. This 
is likely due to the fact that Wilson fermions on a finite lattice introduce a 
mass shift, i. e. an additive mass renormalization, which is not present in the 
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| Wave function renormalization without lattice correction 



Mass function without lattice correction 




Figure 1: The functions Al and Bl, as well as the derived quantities 
Zl = 1/Al and Ml = Bl/Al in lattice units for different lattice settings. 
Statistical errors here and afterwards are smaller than the box size, if not 
visible. 



continuum [56f] . Thus, the continuum renormalization prescription cannot 
be transferred. 

As a consequence, the renormalization constants Za and Zg are inde- 
pendent, and have to be included separately. They are determined using the 
renormalization conditions 



A(fx) 
B(jx) 



1 



m. 



(7) 



with /i = 1.9 GeV. This value is chosen as it reduces discretization er- 
rors, since for all lattice spacings lattice momenta are sufficiently close to it 
available to permit a linear extrapolation of the dressing functions to fi to 
determine the multiplicative renormalization constants Za and Zb- 

Doing so results in the transformation of figure [2] to O where the errors 
of the renormalization constants have been included by error propagation. 



A function with lattice correction 



B function with lattice correction 
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| Wave function renormalization with lattice correction 



Mass function with lattice correction 
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Figure 2: The functions A and B, corrected for leading lattice artifacts, as 
well as the derived quantities Z = 1/A and M = B/A. 



The results show that there are still quite substantial lattice artifacts, which 
will be discussed in turn. For this purpose a bare quark mass is convenient, 
which is neither small nor large with respect to the lattice parameters. For 
this purpose, 1 GeV will be chosen. 

The first lattice artifact investigated is the dependence on the physical 
volume. This is shown in figured] for two different discretizations. Interest- 
ingly, there is almost no volume-dependence on the coarser lattice, while the 
results on a finer lattice show some effect, and in particular a large change 
compared to the coarser lattice. Comparing the case of a mass of 1 GeV 
to the cases of 0.5 and 2 GeV shows that this effect becomes smaller for a 
larger mass (the ratio at small momenta decreases from 1.2 to 1.1) while 
it increases for the smaller mass (to 1.4). This suggests that finite- volume 
effects become strong already at surprisingly large masses, given that mL 
for 2 GeV and the largest volume is about 24. 

Also interesting is the comparison of discretization artifacts at fixed phys- 
ical volume. These are shown in figure [5j It is immediately visible that an 
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| Wave function renormalization with lattice correction 



Mass function with lattice correction 
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Figure 3: The renormalized functions A and B, as well as the derived quan- 
tities Z = 1/A and M = B/A is lattice units for different lattice settings. 
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Figure 4: The renormalized mass function and wave-function renormaliza- 
tion for two fixed lattice spacings and fixed mass and different physical 
volumes. 
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Figure 5: The renormalized mass function and wave- function renormaliza- 
tion for three fixed physical volumes and fixed mass and different discretiza- 
tions. 

improved discretization induces a stronger infrared suppression. However, 
this effect diminishes quickly with increasing physical volume. This effect 
becomes also smaller for larger mass and increases for smaller masses. This 
again emphasizes that for a mass of about 1 GeV, and even for 2 GeV, 
lattice artifacts play a substantial role. 

The wave-function renormalization shows no pronounced dependency on 
either volume or discretization. However, a very close study shows that in a 
statistical significant way the wave-function renormalization becomes flatter 
with improved discretization, see figure [5l while it is almost independent 
of volume. However, one should be wary that in the fundamental case 
different quark actions yield on similar lattice setups quite different wave 
function renormalizations [39l.l40l.l42]. It can therefore not be excluded that 
with a different discretization of the quark propagator the wave-function 
renormalization shows more substantial deviations from one. The impact of 
different quark propagator discretizations on the mass function is, however, 
by far not as severe. 

Summarizing, the lattice artifacts for the lattice settings of table Q] for 
the employed masses are substantial, and induce especially at small mo- 
menta systematic errors exceeding 10%. That is somewhat surprisin g, g iven 
the inherent mass-scale of the quarks, though gluonic correlators [47|] are 
similarly affected. 

Keeping this in mind, the final results are presented in figure [6j 

As already expected from the analysis of the lattice artifacts the wave- 
function renormalization does not deviate substantially from one. In fact, 
it is found that within errors all deviations from one depend on the lattice 
momentum rather than on the physical momenta, and diminish therefore at 
fixed physical momentum with better discretization. Thus, within errors, it 
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I Wave function renormalization with lattice correction 
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Figure 6: The wave function renormalization (left panel) and mass function 
(right panel) for the various masses and discretizations (circles are 0.2 fm, 
squares are 0.15 fm, triangles are 0.1 fm, and upside-down triangles are 0.05 
fm), always from the 24 4 lattices. For the wave-function a constant one is 
also displayed, while for the mass-function fits of type Q are included by 
the full lines. 
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i are dimensionless, al 
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2b 


7 


A 


-<§#)5 
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M(0) 


0.01 


0.2(4) 


in-u 
i - y +1.3 


1 2 +i ' d 


0.12(7) 


0.003/0.01/2xl0 5 


1.6(1.2) 


°- 18 -0.03 


0.1 


0.309(1) 


1.976(5) 


1 2 +ai 


0.35(3) 


9 c+U.3 
z -°-2.3 


i 9+0.2 


1.2(2) 


0.5 
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2.33(3) 


1.64(7) 
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0.560(1) 
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0.7239(8) 
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3.00(9) 


10 


0.361(1) 


0.138(1) 


1.52(2) 


1.968(4) 


1.8(3) 


3.0(2) 


10.9(3) 



appears that the wave- function renormalization is indeed constant. This is 
quite different from the case for fundamental quarks 33j, [35J , and somewhat 
surprising. However, it supports this assumption in several calculations 
using DSEs in the context of technicolor and conformal window studies 

2 EH IS- 

The mass-function shows a significant dependence on the physical mo- 
mentum, and is only rather limited affected by the lattice artifacts. Only for 
the smallest mass of 0.01 GeV unsurprisingly significant lattice artifacts are 
seen, dominating especially at small momenta. It can thus be expected that 
only limited physical information can be obtained from this mass. However, 
even for the largest mass of 10 GeV, very few lattice artifacts are seen, which 
was not expected. 

In all cases the mass function shows an infrared screening mass signifi- 
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cantly larger than the mass implemented at the renormalization point. This 
is usually interpreted as a sign of contributions from chiral symmetry break- 
ing [321 ] . To investigate this further, the mass function were fitted at large 



momenta by two possibilities [57|, [58[ , the one the (up to explicit breaking) 
chiral symmetric case 

(2 \ ~ ' 7 

wln^ + lj , (8) 

where uj and the anomalous mass dimension 7 were used as fit parameters. 
The alternative was the operator-product-expansion motivated form for the 
chirally broken case 

2vr 2 7 -(**) 



M(p) 



2 \i-7' 



p 2 ( 4 In £ 



leaving the characteristic scale A, the chiral condensate (\I n I / ), and 7 as free 
parameters. However, both forms do not satisfactorily fit the curves. This 
is achieved by using a more relaxed form 

,,/ n 2vr2 7 -(#*) 

M(p) = - 



p 2b i In 



2 m A 7 



-7' 



with anomalous exponent 6. This form can be upgraded to a regularized 
version 

m(p) = -^ ; \ -1^ , o) 

with two further parameters a and c, which then fits the quark mass func- 
tion over the whole momentum range. The fit parameters, together with 
the screening mass M(0) are given in table [2j All fits were performed by 
dropping both the lowest and the two largest momentum points on all lat- 
tices, to remove the most dominant lattice artifacts in the fits. Of course, 
7 and (fyfy) in the fit form (|U]) can no longer be strictly interpreted as the 
anomalous dimension and the chiral condensate, except when b ~ 0, which 
is not the case. 

As is visible, the strong finite- volume effects on the smallest mass makes 
the fit completely unreliable, while for the other masses quite stable results 
can be obtained. In particular, the physically most interesting parameters 
are most robust, while the regularization parameters a and c show somewhat 
stronger fluctuations, especially for smaller masses. It is interesting that the 
would-be anomalous mass dimension shows a pronounced dependence on 
the mass, moving towards the unitarity limit of -2 for decreasing mass. 

The results show that the effective screening mass M(0) is indeed about 
1 GeV larger than the bare mass. Interpreting this as a chiral scale, this 
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Figure 7: The vector (left panel) and scalar (right panel) Schwinger functions 
for two fixed lattice spacings and fixed mass and different physical volumes. 
Lines are to guide the eye. 
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Figure 8: The vector (left panel) and scalar (right panel) Schwinger functions 
for three fixed physical volumes and fixed mass and different discretizations. 
Lines are to guide the eye. 

scale is thus much larger than the one for fundamental quarks, which is of 
order 300 MeV [33J . If one takes the value of the chiral condensate seriously, 
which is probably not adequate since the exponent b is far from 1, this would 
suggest rather a chiral scale similar to the one of fundamental quarks. To 
better understand this, it is worthwhile to investigate the quark propagator 
in configuration space, i. e. its Schwinger functions (J2M3D- 

3.2 Results in position space 

Though the momentum-space results do not show substantial lattice arti- 
facts once the most extreme momenta are neglected, this can change in 
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Figure 9: The vector (left panel) and scalar (right panel) Schwinger func- 
tions for the different masses from the largest volume. Lines are the lattice 
Schwinger functions of the corresponding fits ([U]) . 



position space. Therefore, figures [7] and [8] show the analogues of figures 2] 
and [5] in position space, i. e. the dependence on volume and discretization, 
respectively. The severeness of volume artifacts is seen to be stronger the 
better the lattice discretization. While for the coarse lattice this is not too 
severe, it is very substantial for the finer lattice. The discretization effects 
seen in figure [8] are also sizable. However, their size strongly depend on 
the volume, and they become smaller the larger the volume. Thus, the 
most important requirement is a sufficiently large volume. However, for all 
lattice parameters available, sizable effects remain, which will distort any 
quantitative results at least by 10%. 

Keeping these systematic effects in mind, the best choice appears to use 
the largest physical volume. The behavior for the different masses is then 
shown in figure The first observation is the presence of a dip in all cases, 
except for the largest mass. This dip in the plotted absolute value of the 
Schwinger functions originate from a sign change in the Schwin ger function 
itself, thus implying a positivity violating spectral function 43|, |43]. Hence, 
quenched adjoint quarks are not part of the physical spectrum, similar to 
gluons [43]. The position of this dip depends strongly on the mass, and it 
is generally at shorter distances for the vector Schwinger function than for 
the scalar Schwinger function. 

The corresponding (lattice) Schwinger functions of the fits ([9]) describe 
the behavior rather well, except for the lightest mass. Interestingly, they 
signal the presence of more zero crossings, though this cannot be unambigu- 
ously confirmed by the lattice data itself due to the large statistical fluctua- 
tions. Such an oscillatory behavior would indicate conjugate complex poles 
as the analytical structure 43l. 1471] . This would be in-line with the commonly 
assumed structure for the fundamental quark propagator 32|, |33|] . 
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4 Unquenched 

Before discussing the results for the unquenched adjoint quark propagator 
some particularities of the unquenched configurations have to be addressed. 
The first one is with regard to the scale. There is no (known) equiva- 
lent to the unquenched theory in nature. The results on the bound state 
spectrum 18H21] furthermore suggest that the physics is completely differ- 



ent from the QCD case, with the lightest glueball always being lighter than 
the lightest hadron. Especially, breaking of chiral symmetry in the chiral 
limit could not be established, and the theory is expected to be walking or, 
possibly, quasi-conformal or even conformal in the chiral limit [1J, I la . 121 ] . 



This makes the scale setting quite arbitrary. Based on technicolor phe- 



nomenology, in the previous work 30(] the lightest state in the channel ++ , 
the dilaton, being here the scalar glueball [19j], has been used in to set the 
scale, assuming a mass of 2 TeV. This scale provided a rather smooth be- 
havior of gluonic correlators with the quark mass, and appears therefore 
reasonable. It also coincides, within errors, with the scale determined using 
the conventional fundamental string tension. It is this approach which will 
be followed her to set a physical scale. In addition, to avoid any bias in- 
troduced by this scale setting, all final results will also be shown in lattice 
units. 

However, the value of 2 TeV appears questionable. Based on the consid- 
erations in [53|, [54( , it should be expected that the scalar glueball/dilaton 



state actually mixes with the observed [5l|, l52| Higgs-like particle, if it has 
not even to be identified with it. Based on these considerations, this state 
should therefore be attributed rather a mass of 125 GeV, which will be done 
here. 

Finally, there is the question of the renormalization scheme. Without 
explicit referral to the quarks, previously it sufficed to just use the bare quark 
mass [30j . However, for the unquenched case most of the bare masses are 
actually negative 18h21|. especially towards the chiral limit, and therefore 



unsuitable for the renormalization prescription. To address this problem, 
the PCAC quark mass |la |. being a positive quantity, will be used here 
instead of the bare mass for the renormalization condition ([7]). Since this 
mass is only available over a limited range, it was extrapolated/interpolated 
linearly for all other cases, particularly towards large bare quark masses. 
This mass was also used for the lattice correction Q. Since in the limit 
a — > this choice drops out in the renormalization procedure, any deviation 
from an optimal choice to remove la! ttice artifacts only surfaces again as a 
lattice correction. 

Note that the configurations available to us are only at one fixed bare 
coupling /3 = 2.25, though with varying bare quark mass, and only for a 
very limited range of volumes. This limits the amount of systematic analysis 
possible severely compared to the quenched case. 
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Figure 10: The vector (top left panel) and scalar (top right panel) Schwinger 
functions and the wave function renormalization (bottom left panel) and 
mass function (bottom right panel) for different volumes for am = —1.175, 
i. e. a PCAC mass of 0.0678. 



4.1 Impact of the wrong phase 

A serious problem surfacing is that a spatial center-breaking transition is 
observed [19| when decreasing the bare lattice mass towards the chiral limit 
of am w — 1.2 [la ]. It is observed for a spatial extent of N s = 8 at bare quark 
masses below am = —0.975, for N s = 12 below am = —1.05, and for N s = 16 
to set in between am = —1.05 and am = —1.15 18j. The spatial dependence 
suggests that changing the bare quark mass at fixed gauge coupling (3 also 
changes the lattice spacing, and especially reducing it toward the chiral limit. 
This is supported by other observations in the gluonic sector [30J ■ Thus the 
physical volume shrinks when lowering the bare quark mass at fixed (3, and 
at a certain physical extent a phase transition occurs, just like in Yang-Mills 
theory. Therefore, in the following this transition will be assumed to be a 
finite- volume artifact. 



This did neither influence significantly gluonic correlators 30J], nor did 
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it distort the physical spectrum in lattice units lallS)]- In contrast, and 
thus surprisingly, we find here that it does alter the quark dressing func- 
tions, especially the B function, qualitatively. Especially, it induces a sign 
change at a finite momentum. This is illustrated in figure [TUJ While the 
vector Schwinger function and the wave-function renormalization appear 
still mostly normal, this is dramatically different for the scalar part. The 
scalar Schwinger function increases at first. This is not yet a problem, as 
it does so also in the quenched case after becoming negative, see figure El 
The consequence for the mass function is much more severe: It has a sign 
change. This originates from a sign change of the uncorrected B function. 
This behavior first increases with volume from N s = 8 (where the function 
is st! ill positive) to N s = 12. Reducing the aspect ratio from 1 : 2 to 1 : 1, 
this effect is counteracted. It is also reduced when going to N s = 16. The 
latter is expected for a finite- volume artifact. Unfortunately, there a not 
enough different volumes available to check this explicitly close to the phase 
transition. Therefore, it can only remain a conjecture that this is a pure 
finite- volume artifact. 

When instead of the bare quark mass a sufficiently larger quark mass 
is used when determining the quark propagator, i. e. a larger bare valence 
quark mass than sea bare quark mass, the B function becomes once more 
positive. However, since it is not clear what kind of physics such a partially 
quenched theory actually describes, this will not be explored further here. 

This serious lattice artifact therefore prevents any physical interpretation 
of the quark propagator inside the wrong phase, and limits the following 
investigation to the cases with volumes large enough to have at least the 
possibility to be in the infinite-volume phase. 

4.2 Results in momentum space 

The results for the wave function renormalization and mass function in both 
lattice and physical units in the relevant bare mass range are shown in figure 
[TT1 Again, a fit with the form ([9|) is possible, and the parameters are shown 
in table [3l However, since the exponent b is close to zero for small quark 
masses, and 7 is always close to one, results for the constants a, c, and A 
are quite unreliable and very sensitive to statistical fluctuations. They are 
therefore not listed, since their usefulness is very limited. 

The wave-function renormalization is, as in the quenched case, for all 
practical purposes indistinguishable from one. For the mass function, how- 
ever, there are quite a number of differences compared to the quenched case. 

The first is that the mass function is much slower decreasing towards 
the ultraviolet. Especially, for smaller quark masses the exponent b is essen- 
tially zero, while the anomalous exponent 7 is close to one. In fact, for the 
smallest quark masses accessible without phase transition the mass function 
is essentially constant. Also, the effective mass M(0) is much less enhanced 
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Figure 11: The wave function renormalization (left panels) and mass func- 
tion (right panels) in lattice units (top panels) and physical units (bottom 
panels) for different masses and volumes. For the case of lattice units also 
fits of type ([9]) are shown. Renormalization is performed at o/i = 1. 

compared to the mass function at the renormalization point. This ratio is, 
within errors, even one for the smallest quark masses. This indicates that for 
small quark masses there is no spontaneous chiral symmetry breaking left, 
and only explicit breaking remains. These statement holds true independent 
of the units. 

In physical units, in contrast to the case of the gluonic correlation func- 
tions [30(, the mass function does not collapse onto a universal curve. Given 
that the amount of chiral symmetry breaking is reduced with decreasing 
mass, and therefore also the effective mass M(0) decreases quicker, this is 
not surprising. Interestingly, the chiral condensate in physical units stays 
more or less constant, as do the PCAC masses in physical units, as soon 
as the PCAC mass reaches a certain level. At this level, the curves again 
collapse to a certain degree to the same curve. 

These results are not easily interpreted. It appears that moving to 
smaller masses reduces the amount of spontaneous chiral symmetry break- 
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ing, until it essentially vanishes. This would be in line with the expectation 
that the theory develops a (quasi- )conformal infrared fixed point in the chiral 
limit jjj&lSi. 

However, whether this is indeed an effect of reaching the chiral limit, or 
whether this is a finite- volume artifact due to the vicinity of the spatial phase 
transition is less easy to understand. Since the effect seems to be rather 
volume insensitive at fixed quark mass, this may be possibly assigned to the 
former option. However, without much more systematic investigations, this 
could also be wishful thinking. Thus, this result is more an encouragement 
to extend the systematics, rather than a final result. 

With the exponent b being close to zero, the form of the propagator is 
close to the form ([8]). This would permit to interpret 7 at small quark masses 
indeed as the anomalous dimension, which would then be unexpectedly large, 
and close to one. This is in stark contrast to indirect determinations possible 
59j, l60|] from gauge-invariant observables, like the hadronic spectrum [18] or 
the spectrum of the Dirac operator 31> 60]. Whether this is a lattice artifact, 
a (residual) scheme dependency given that Za 7^ Zb due to lattice artifacts, 
or has another reason is a question which has to be explored. 

However, there is, of course, a certain arbitrariness in the fit form (|9|). 
Given that the largest accessible momenta are not too large compared to all 
other scales, the regularization can influence the result. If, e. g., a fit form 
generalizing (|HJ) of type 



M(p) 



Uln j2 +k 



-lv 



is used instead, the obtained anomalous exponent 7 P is found to be of order 
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Figure 12: The vector (left panels) and scalar (right panels) Schwinger func- 
tions in lattice units (top panels) and physical units (bottom panels) for 
different masses and volumes. For the case of lattice units also fits of type 
([9]) are shown. Renormalization is performed at a\x = 1. The symbols are 
the same as in figure [TT1 

0.2-0.3 for small quark masses [6l(. This is in much better agreement to 
other available determinations, e. g. [18|, |3l|, l60f] . Without chiral symmetry 
breaking, this is also the more adequate form [571. l58l]. 

Of course, given that b ~ in ([9]) for small quark masses, it is possible 
to identify the exponents of the logarithms, j p = 1 — 7. This would yield 
then a j p close to zero, and therefore smaller than other determinations. 
Given this relevance of the fit form, this should be regarded as a systematic 
error. Hence, the result is, given all the assumptions, that 7 ~ 1 if chiral 
symmetry is spontaneously broken in the chiral limit, and 7 < 0.3 if not. 



4.3 Results in position space 

The Schwinger functions, shown in figure fT2| show no significant changes, 
compared to the quenched case. The bending in all cases shows that the 
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dynamical adjoint quarks do not posses a positive-definite spectral function, 
and therefore do not belong to the physical state space. At larger masses, 
the zero crossing is still explicitly observable, but this is no longer possible 
at small quark masses. Also, in both lattice and physical units, the position 
of the first zero crossing moves towards larger times the smaller the mass 



Though the violation of positivity is not sufficient for confinement 47( , 
it implies that the quarks are not observable particles. Thus, they must be 
either bound in bound states or confined, for all the masses investigated 
here. Furthermore, the behavior is neither a (approximate) power-law nor 
close to one, as would be expected if the theory would be ( quasi- )conformal 
[621 ]. In fact, the fits of type © describe the Schwinger functions quite well. 

Combining both results, speculatively the behavior of the adjoint quarks 
in the limit of zero quark mass is best characterized as that of confined 
quarks without spontaneous chiral symmetry breaking. This underlines that 
adjoint quarks are substantially different from fundamental quarks. 

5 Summary 

Summarizing, the investigation of the adjoint quark propagator presented 
provided a number of interesting observations, though further systematic 
studies are necessary for any final statements. In any case, adjoint quarks 
are quite different from fundamental quarks. Especially their wave-function 
renormalization appears to be essentially identical to one, in strong contrast 
to the fundamental case 33|, [35(. Furthermore, in the chiral limit the results 



suggest the possibility that spontaneous chiral symmetry breaking vanishes 
in the chiral limit, at least for 2 colors and 2 flavors. At the same time 
the Schwinger functions show clearly a violation of positivity. Thus adjoint 
quarks do not belong to the physical spectrum, irrespective of whether in 
the quenched case or not. Finally, the mass function can be rather well 
described by the form ([9]). The interpretation in terms of an anomalous 
dimensions hinges crucially on the presence of spontan! eous chiral symme- 
try breaking in the chiral limit. If this breaking is absent, the anomalous 
dimension is small, in agreement with other determinations, otherwise it is 
large, close to one. Speculatively, this and other indications presented here, 
could be taken as a sign for the absence of chiral symmetry breaking in the 



chiral limit, in line with [14l43ll|. However, a much larger range of momenta 
would be necessary to be sure to have eliminated all residual systematic 
uncertainties introduced in the fitting process. Thus again a significantly 
enlarged investigation is called for. 

In total, the results clearly show that adjoint quarks and fundamental 
quarks are very different, already for the quenched case, and become even 
more different towards the chiral limit. But much more systematic investiga- 
tions will be needed to provide the necessary systematic reliability to make 
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this a final statement. To really understand the chiral limit will probably 
require a chiral extrapolation of the dynamical case, which can be achieved 
using, e. g., functional methods, like done already for fundamental quarks 
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